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VARIATIONAL  STUDY  OF  NONLINEAR  SPLINE  CURVES 


by 

K.  H.  Lee  and  G.  E.  Forsythe 
Stanford  University 

1.  Linear  and  nonlinear  splines.  Let  A  be  a  finite  ordered 
set  of  points  in  the  euclidean  plane,  with  cartesian  coordinates  (x^,yr)  , 
r  1,  ...,n-l  ,  through  which  it  is  desired  to  pass  a  smooth  curve.  An 
old  technique  in  drafting  is  to  use  a  mechanical  spline  to  form  a  smooth 
curve  C  that  contains  A  .  In  the  present  day  of  automatic  plotters, 
numerically  controlled  milling  machines,  and  so  on,  it  is  more  important 
to  find  a  mathematical  or  computational  representation  of  a  suitable  C 
than  to  draw  it.  Thus  one  uses  some  mathematical  model  of  the  mechanical 
spline. 

By  far  the  most  widely  used  model  is  a  linear  (cubic)  spline,  suitable 

when  the  curve  C  in  some  x-y  coordinate  system  is  the  graph  of  a 

function  f  ,  so  that  y  =  f(x)  ,  xQ  <  x  <  xfi  .  Assume  that 

xQ  <  x1  <  . . .  <  xn  .  The  linear  spline  can  be  defined  as  the  unique 

function  f  for  which 
x 

f  n  r(x)?dx  (i) 

xo 

is  minimized  among  all  twice  continuously  differentiable  functions  assuming 
the  value  yr  at  xr  (r  -  1, 2,  ...,n-l)  .  (According  to  this  definition, 
f  will  satisfy  the  natural  end  condit ions  f"(x0)  =  f"(xn)  --  0  .  There 
are  alternative  treatments  of  the  end  conditions.) 


The  linear  spline  so  defined  turns  out  to  be  a  (usually)  different 
cubic  polynomial  in  each  interval  (x^  ^,x^)  *  with  matching  values, 
derivatives,  and  second  derivatives  (and  hence  curvatures)  at  each 
interior  node  x^  (r  =  1, 2, ...,n-l)  .  The  spline  will  actually  be  a 
straight  line  seonent  for  x,.  <  x  <  x,  and  x  n  <  x  <  x 

The  theory  of  linear  splines  has  grown  enoimously  in  the  last  decade, 
and  these  curves  and  various  linear  generalizations  have  both  practical  and 
theoretical  importance  in  the  approximation  of  known  functions,  solutions 
of  differential  equations,  and  so  on.  The  reader  can  find  an  exposition, 
with  generalizations  and  applications,  in  Ahloerg,  Nilson,  and  Walsh  (1]. 

Linear  splines  are  invariant  under  linear  changes  in  the  y-coordinate 
alone,  as  Podolsky  and  Denman  [9]  point  out.  Thus  they  are  suited  to 
such  problems  as  the  interpolation  of  data,  where  x  and  y  have 
different  meanings.  On  the  other  hand,  linear  splines  are  not  invariant 
under  rotations  of  the  x-y  coordinate  system,  and  hence  are  not  well 
suited  to  the  interpolation  of  geometrical  points  in  the  euclidean  plane. 
Moreover,  linear  splines  cannot  be  used  directly  to  define  a  closed  curve  ‘ 

C  in  the  x-y  plane. 

For  the  purposes  of  interpolating  points  in  the  euclidean  plane  it 
is  appropriate  to  find  a  mathematical  model  which  is  invariant  under  all 
similarity  transformations.  The  model  we  treat  is  sometimes  called  the 
elastica,  but  we  shall  refer  to  it  as  a  nonlinear  splin  ■.  As  a  preliminary 
to  work  on  actually  computing  nonlinear  splines,  we  have  investigated  their 
precise  definition,  including  variational  properties,  defining  equations, 
continuity  conditions,  and  end  conditions,  both  for  open  and  closed  curves. 
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The  term  nonlinear  spline  is  used  variously  in  the  literature: 

(a)  If  the  integrand  f'*(x)  of  (1)  is  multiplied  by  a  nonconstant 
weight  factor,  sometimes  the  function  that  minimizes  the  altered  problem 
is  called  a  nonlinear  spline,  (b)  Suppose  one  is  given  a  function  <p(x) 
to  be  approximated  by  a  linear  cubic  spline  passing  through  n-1  points 
(xr,«p(xr))  ,  and  that  the  n-1  abscissas  xr  are  varied  until  the  spline 
best  approximates  <p  in  some  given  norm.  Sometimes  the  result  is  called 
a  nonlinear  spline.  However,  in  both  (a)  and  (b)  above  the  splines 
satisfy  a  linear  differential  equation  in  each  interval,  whereas  our 
nonlinear  splines  satisfy  a  nonlinear  differential  equation. 

We  do  not  claim  that  computing  nonlinear  splines  will  necessarily  be 
an  economical  way  to  interpolate  points  in  the  x-y  plane.  Moreover, 
nonlinear  splines  are  not  invariant  under  linear  changes  in  the  y-coordinate 
alone,  so  that  they  seem  ill-suited  to  the  interpolation  of  data  where  x 
and  y  are  unrelated  . 

We  have  been  interested  only  in  studying  as  carefully  as  we  could  the 
mathematical  nacure  of  these  nonlinear  splines.  In  this  paper  we  present 
a  variational  treatment  of  nonlinear  splines,  emphasizing  the  natural 
boundary  conditions  of  the  problem.  We  believe  that  our  treatment  of 
the  closed  nonlinear  spline  may  be  new. 

2.  Previous  work.  In  the  theory  of  elasticity,  our  mathematical 
model  of  the  mechanical  spline  is  called  a  thin  beam  or  elastica,  and  its 
treatment  dates  back  to  James  and  Daniel  Bernoulli,  Euler,  Kirchhoff, 
and  others.  The  history  and  theory  are  summarized  by  Love  [7].  None  of 
t  iose  treatments  dealt  directly  with  the  use  of  the  nonlinear  spline  to 
interpolate  points,  and  there  was  little  discussion  of  closed  splines. 
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The  earliest  discussion  that  we  have  seen  in  print  of  the  use  of 
nonlinear  splines  for  interpolation  is  that  of  Hirkhoff  and  do  boor  |2|. 
That  paper  refers  to  excellent  laboratory  reports  by  Fowler  and  Wilson  [t] 
and  by  Birkhoff,  Burchard,  and  Thomas  [5],  Glass  [5]  briefly  describes 
computations  of  open  nonlinear  splines  in  cartesian  coordinates. 

Hosaka  [6]  describes  the  generation  of  nonlinear  splines  on  a  digital 
differential  analyser.  Woodford  [12]  describes  an  iterative  procedure 
for  interpolation  with  open  nonlinear  splines  that  is  much  faster  than 
Glass's  algorithm;  he  also  works  with  cartesian  coordinates. 

In  his  Ph.D.  dissertation.  Medium  [8]  discusses  the  nature  of 
nonlinear  open  splines,  again  using  a  cartesian  coordinate  system  rotated 
to  a  convenient  local  orientation.  He  also  gives  an  algorithm  for 
computing  an  approximation  to  the  nonlinear  spline  by  a  succession  of 
circular  arcs  meeting  with  a  continuous  tangent  but  discontinuous 
curvature . 

3.  Basic  concepts  of  bending  theory  of  thin  beams.  Of  all  the 
curves  that  pass  in  turn  through  the  ordered  set  A  of  points  (,x^,yr) 
mentioned  in  the  introduction,  we  shall  consider  as  admissible  only  those 
whose  tangent  direction  is  continuous  everywhere,  and  whose  curvature  is 
piecewise  continuous,  with  discontinuities  in  curvature  permitted  on 
any  finite  set  of  points.  A  plausible  suggestion  for  the  smoothest  of 
these  admissible  curves  is  that  the  integral  of  the  square  of  :he  curvature 
with  respect  to  arc  length  should  attain  a  minimum.  This  comprises  a 
simple  representation  of  the  concept  of  a  curve  passing  through  the  points 
with  minimum  total  bend  amplitude,  and  will  be  utilized  in  the  form  of  the 
necessary  condition 
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2 

K  ds  =  0 


(2) 


f 

n 


where  >v  is  the  curvature,  s  is  arc  length,  *n"*o  **  the  total  length 

of  the  curve,  and  o  is  the  symbol  for  variation.  The  integral  in  (2)  is 
proportional  to  the  strain  energy  in  a  bent  spline  according  to  Euler- 
Bernoulli  beam  theory,  and  we  show  in  Section  1+  that  (2)  is  the  variational 
form  of  the  conditions  of  equilibrium  for  the  spline  with  forces  applied 
only  at  the  support  points.  It  seems,  therefore,  appropriate  to  investigate 
spline  interpolation  in  terms  of  mechanical  bending  theory,  and  it  will  be 
shown  in  the  present  paper  that  this  approach  does  lead  to  the  introduction 
of  variables  which  are  particularly  convenient  for  interpreting  spline 
interpolation,  and  perhaps  also  for  computing  splines. 

Bernoulli -Euler  theory,  as  described  in  detail  by  Love  [7],  is  the 
simplest  form  of  beam  theory,  and  considers  only  bending  deformations, 
neglecting  shear  deformations  and  stretching  of  the  center  line  of  the 
beam.  Guch  an  approximation  is  satisfactory  for  Deams  with  cross-section 
dimensions  small  compared  to  the  span  between  supports,  as  clearly  applies 
for  splines.  Such  restricted  deformations  are  introduced  by  requiring 
that  plane  sections  normal  to  the  center  line  in  the  undeformed  state 
remain  plane  and  normal  to  the  deformed  center  line,  and  that  the  center 
line  of  the  beam  does  not  stretch. 

The  forces  and  moments  on  a  beam  element  are  shown  in  Figure  1, 
where  M  is  the  bending  moment,  S  is  the  shearing  force,  and  P  is 
the  longitudinal  force.  (The  convention  is  P  >  0  for  tension,  and 
P  <  ">  for  compression.)  The  assumptions  about  deformation  mentioned 
above,  combined  with  Hooke's  law  relating  stress  and  strain,  yield 


M  =  E  I  K  ,  (3) 

where  E  is  Young's  modulus  of  elasticity,  I  is  the  second  moment 
of  the  section  about  the  axis  of  bending, 


K 


&e 

ds 


> 


(M 


and  6  is  the  angle  between  the  beam  and  the  x-axis.  The  equations  of 
equilibrium  for  each  unloaded  span  between  supports,  which  are  deduced  in 
Section  k,  are  as  follows: 


for  moments: 


=  o  ; 


for  normal  forces: 


dS 

dF +  P*  =  0  5 


for  longitudinal  forces: 


dP 

d^-S«  =  0  * 


It  is  convenient  to  work  in  terms  of  reduced  force  variables: 


M  = 


El 


=  K  ,  S  = 


JL  p  _  _P 

EI  '  El 


Then  (5)  and  (6)  give: 


5  d  K 

S  =  ’  di  ’ 


P  =  — 


2 

1  d  k 


K  .  2  ’ 

ds 


and  (7)  becomes 


(5) 


(6) 


(7) 


(8) 


(9) 

(10) 
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0 
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(11) 


_d_  r  _L_  djc 

ds  1  K  .  2  2 

ds 


or 


1 

K 


d'V 


ds 


K 

n 


C 


r-l  1 


(12) 


where  cr_2  i-s  a  constant  of  integration  for  the  r-th  span. 

We  must  now  consider  the  boundary  conditions  at  the  supports  which 

constrain  the  spline  to  pass  through  the  required  points  Q  of  A  . 

Tne  least  constraining  such  support  is  a  freely  rotating  sleeve  attached 

i/O  the  point  that  permits  free  rotation  of  the  spline  and  free 

sliding  through  the  sleeve.  The  only  support  force  is  therefore  normal 

to  the  sleeve,  and  this  does  no  work  on  a  possible  motion  of  the  spline 

tiirough  the  sleeve.  A  more  constraining  support  would  be  a  pin  through 

the  spline  which  permits  free  rotation  but  no  sliding,  or  a  pin  with 

rotation  prevented.  In  none  of  these  cases  is  work  done  by  support  forces, 

since  either  a  force  (or  moment)  component  is  zero,  or  the  associated 

motion  is  zero,  and  such  supports  are  termed  workless  constraints. 

Figure  2  shows  a  spline  passing  through  freely  rotating,  sliding 

sleeves  at  Q  ,  ...,Q  ,  ,  where  Q_  and  Q,  are  the  free  ends  of  the 

1  n-1  o  ^n 

spline.  The  configuration  of  the  spline  could  be  analyzed  using  the 

equations  given  above,  but  a  simpler  and  more  revealing  approach  for  our 

purposes  is  to  observe  that  this  spline  forms  a  conservative  mechanical 

system  with  potential  enemy  given  by  the  strain  energy  of  the  spline, 

i 

,  n  o 

u  =  ^i-r/.Tjds  ;  (13) 

i 

{ t 
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there  are  no  other  contributions  to  U  ,  since  the  external  forces  are 
all  worklesr.  The  theory  of  conservative  systems  [11]  tells  us  that  at 
a  stable  equilibrium  configuration  of  the  spline,  the  energy  (15)  is  a 
local  minimum,  which  implies  (2)  for  a  uniform  spline  with  El  constant. 
Moreover,  any  constraint  added  to  the  system,  such  as  changing  a  freely 
sliding  sleeve  to  a  pin  support  that  prevents  sliding,  will  either  increase 
tlie  potential  energy  in  the  corresponding  equilibrium  configuration,  or 
Leave  it  unchanged  if  the  added  constraint  happens  to  be  compatible  with 
the  configuration.  Thus 
l 

n  o 

J  <  ds  (lb) 

*0 

will  also  exhibit  a  local  minimum  in  the  configuration  shown  in  Figure  2, 
relative  to  variations  of  the  constraints.  Note  that  the  free  ends 
and  ,  with  no  forces  or  moments  applied,  also  provide  workless  boundary 
conditions,  and  any  constraint  on  their  freedom  of  motion  will  increase 
the  energy  expression  (lb).  Thus  a  local  minimum  of  the  integral  (lb) 
corresponds  to  free  ends  and  freely  rotating  sliding  constraints  at 
.  ..,Q  ^  .  This  cannot  bo  a  global  minimum  in  the  space  of  all 
configurations,  since  a  lower  value  of  the  integral  in  (3b)  can  be  achieved, 
as  pointed  out  in  [2],  by  introducing  large  loops  between  supports,  which, 
of  course,  modify  the  topology. 

In  the  next  two  sections  we  deduce  the  least  constraining  support 
conditions  for  the  spline  passing  through  the  points  Q  , ...,Q  ^  by 

seeking  the  minimum  of  the  integral  (lb)  directly  through  analysis  of  the 
variational  problem  (2),  and  deduce  the  natural  boundary  conditions  that 


yield  this  minimum  energy  configuration.  Although  this  approach  simply 


T 


reproduces  the  minimum  constraint  conditions  shown  in  Figure  2,  and 
anticipated  above  on  the  basis  of  cons  ^rvative  system  theory,  it  is 
independently  useful,  since  it  permits  investigation  of  the  closed  spline 
problem  in  Section  (<.  The  latter  problem  cannot  be  treated  directly  by 
the  theory  of  a  constrained  conservative  system,  because  we  must  consider 
the  effect  of  variable  arc  length  for  the  closed  curve,  and  this  changes 
the  system  more  than  simply  by  imposing  a  constraint. 

1*.  Deductions  from  the  variational  statement.  We  consider  the 
variational  statement  (2)  with  integration  limits  and  l  for  a 

curve  constrained  to  pass  through  the  points  Q,^, ... ,Qn  ^  with  end 
points  Q,0  and  Q,  .  The  points  ...,Qn  correspond  to  values  lQ,  ...,1^ 
of  the  arc  length  s  .  In  the  present  section  we  do  not  consider  end 
conditions.  Thus  we  do  not  care  whether  the  curve  is  open  (as  in 
Figure  or  closed  (Qq  and  coincide).  We  shall  prove  in  this 

section  that  if  such  a  curve  satisfies  (2)  --  and  is  hence  a  spline  in 
our  sense  --  then  the  spline  is  the  position  of  a  thin  beam  satisfying 
equations  (5)  -  (7)  of  Section  % 

Because  of  the  constraints,  (2)  takes  the  form 
l 

&  ,  r  2 

°  21  j  K  ds  =  0  (15) 

r=1  t 

r-1 

The  fact  that  the  spline  passes  through  the  points  ...,Q  prescribes 
the  following  constraint  conditions  for  r  -  1, ...,n  : 


10 


0 


> 


I’  cos  8  ds  -  x  +  x  .. 
J  r  r-1 


r-1 

I 


J  sin  fi  ds  -  yr  +  yr_1  =  0  , 


(J'«) 


(1/b) 


r-1 


where  Q  has  the  coordinates  (x  ,y  )  .  Note  in  (1/a),  (l6b)  that  x 


r  r' 


and  are  prescribed  numbers  for  r  =  1,  ...,n-l  ,  whereas  ,  yQ  ,  x 


and  y^  are  free  to  vary. 


We  follow  the  standard  techniques  of  the  calculus  of  variations  and 
introduce  Lagrange  multipliers  Xr  and  ^  for  (1/a)  and  (1/b), 
respectively  (r  =  1, ...,n)  .  We  take  care  of  the  constraints  (1/a),  (1/b) 


by  seeking  a  stationary  value  of  the  functional 


n  r  „  n 

l  J-  *  t  V-1 


r=l 


r=l 


r-1 


J 


cos  wds  - 


x  +  x  , 
r  r-1 


’r-1 


n 

+  £  M.r.x 

r^l 


J  sin  0ds  -  yr  +  y^ 


r-1 


n 


r=l 


*■  ^ 

j  (x  +  \r-1  COS  6  +  ur_1  sin  e)d£ 


r-1 


+  X  . (x  -  x  )  +  j.  (y  -  v  ) 
r-1  r-1  r  *r-lv  r-1  y' 


(17) 
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with  respect  to  a  general  smooth  variation  60(s)  ,  and  variations 


,  by^  ,  bx  ,  fry  ,  combined  with  sliding  through  the  pivots 


6s(Q  )  =  &f 
'  r'  r 


( r  -  1,  •  •  .,n-l)  • 


Setting  the  variation  of  (17)  to  zero  and  integrating  by  parts,-'  we  get 
the  form 


f  It 

n  I  r  r 

z  (  m  * ; 

r'1  I  t*  , 

r-l  r-l 


J*  ’  Vl  Sin  6  +  ‘‘r-l  C0S  e]be  ds 


+  [ +  \  i  cos  6  +  a  ..  sin  9]  52 

r-l  ^r-1  ^  -  r 

r 


-  +  Xr.x  cos  0  +  tir_1  sin  fl]  +  6fr_i 

*r-l 


*  \>  %  *  MO**-,  -  Vl  -  ^n.i  &yn  -  0 


-  In  integrating  by  parts,  we  assume  that  the  curvature  k( s)  of  the 
minimizing  curve  is  continuously  differentiable  in  each  interval 
t  _  ^  <  z  <  <r  .  If  the  curvature  <(c)  of  the  minimizing  function 

is  assumed  only  to  be  piecewise  continuous,  but  0(s)  is  continuous, 
then  it  can  be  proved  by  a  different  argument  based  on  a  lemma  of 
lu  Bois-Reynond  that  -(s)  is  in  fact  continuously  differentiable 
in  each  interval.  This  justifies  our  introduction  of  the  broad  class 
of  admissible  curves  a*  the  start  of  Section  3. 


/(>  ,  /  in  tne  following  we  mean  the  limiting  values  2  +0 


an  i  2  - 


The  integral  term  of  (19)  yields  for  r  =  1, ...,n  : 


-  2  —  -  \  ,  sin  0  +  u,  ,  cos  v  =  u  ,  1  -  ^  s  ^  1  ,  li.'ui 

ds  r-1  ^r-1  r-1  r  ’  v  ’ 

which  can  be  integrated,  using  equations  analogous  to  (l^a)  and  (li  b) 
for  an  open  interval,  giving 

k(s)  =  *(1*^)  -  (y-yr_1)  +  (x-xr_1)  ,  (r  =  l,  (21) 

Identifying  k  with  M  ,  as  in  (8),  we  see  that  (21)  comprises  a  moment 
relation  for  the  part  of  the  spline  between  the  arc  lengths  and  s  , 

as  illustrated  in  Figure  Thus  the  Lagrange  multiplier  factors 
\r  -j/2  and  ^  /2  are  simply  the  force  components  acting  on  the  spline 

at  f*  ^  (r  -  1,  ...,n)  .  By  equilibrium  considerations,  these  same  force 
components  can  be  considered  to  act  on  any  section  of  the  spline  with 
f <  s  <  ,  so  that,  taking  components  along  and  normal  to  the  spline, 

the  tensile  force  P  and  shear  force  S  are  given  for  ^  '  s  <  l ^  by 


cos  0=0  ,  /  1  <  s  <  I  , 


r-1  ^r-1 

cos  0  -  sin  0  , 


"r-1  .  ^r-l  . 

—  sin  0  -  —77—  cos  0  . 


(22a) 


(22b) 


Differential  equation  (20)  can  be  alternatively  integrated  by 


writing 


dK  d  K  d0 
ds  _  d0  ds 


de  ' 


whence,  in  view  of  (22a), 


rrr  +  P  =  c 


(t  -j  < 
r-1 


where  cf  ^  is  an  integration  constant. 


1; 


Note  that  (PO)  and  (P2b)  yield  (9).  Differentiating  (20)  with 
respect,  to  s  and  using  (22a)  give  (10).  Finally,  (12)  and  (11) 
follow  from  (2j5) .  The  basic  equations  (5),  (6),  and  (7)  simply  express 
(0)  -  (12)  in  different  variables,  and  hence  the  equilibrium  equations 
(5)  -  (7)  are  consequences  of  the  variational  statement  (2). 

It  could  conversely  be  proved  that  the  satisfaction  of  equations 
(;>),  (•■),  and  (7)  implies  that  the  variational  condition  (2)  holds.  Thus 
the  variational  condition  (2)  and  the  thin  beam  equations  (5)  -  (7)  provide 
equivalent  foundations  for  the  theory  of  nonlinear  splines. 

Since  be  is  a  continuous  variation, 

be(f-)  =  &0(l+)  (r  =  l,...,n-l)  ,  (21*) 

and  the  first  term  of  (19)  then  demands  that 

K(*‘)=  K(t*r)  (r  -  l,...,n-l)  .  (25) 

In  view  of  (18),  (22a),  and  the  terms  in  (19)  containing  5 , 
we  then  find  that 

P(*y)=P(*r)  (r  =  1, . .  .,n-l)  .  (26) 

5.  The  open  spline.  For  the  configuration  shown  in  Figure  2  with 
free  ends,  6x^  ,  byQ  ,  60(1^)  ,  bx^  ,  6y^  ,  &0(*n)  are  arbitrary 
variations.  Hence  the  first  iud  last  terms  in  (19)  demand  that 

K(l0)  =  =  \  =  ^0  =  Xn-1  =  ^n-1  =  0  *  ^ 2 ^ 

Thus  from  (22a,  22b)  the  end  conditions  become 

,(*0)  P(V  =S(/0)  «*(£„)  -P(ln)  =S(!n)  =0  . 

lb 


(28) 


Thus  the  variational  condition  (15)  implies  that  the  open  spline 

satisfies  the  natural  boundary  conditions  (25),  (26),  (?8),  which 

are  precisely  the  conditions  associated  with  the  least  constraining  supports 

depicted  in  Figure  2  and  discussed  in  Section  3.  In  view  of  'hose 

relations,  (23)  holds  for  the  entire  spline  <  s  <  i  with 

c  ,  =  0  for  all  r  : 
r-1 


%•  +  P  =  0 


and  hence  the  differential  equation 


+  i.  ,  o 
,22  U 


is  valid  as  a  special  case  of  (12)  for  the  threaded  spline  with  free 
ends.  This  equation  has  been  given  by  Birkhoff  et.  al.  in  13 ]•  Note 
that,  in  view  of  (25),  (30)  requires  d  k/ds *  to  be  continuous  across 
supports,  although  in  general  dk/ds  is  discontinuous,  because  the 
lateral  support  force  changes  the  shear  force  S  ,  which  satisfies  (9)  . 

We  wish  to  emphasize  that  our  equations  apply  to  any  spline  curve 
that  satisfies  the  constraints  of  the  problem,  no  matter  what  its 
topology.  As  is  pointed  out  in  [2],  there  may  be  sets  of  nodes  A  for 
which  no  spline  exists  and,  if  any  spline  exists  for  A  ,  there  may  exist 
others  satisfying  the  same  constraints,  with  different  numbers  of  loops 
between  some  adjacent  pair  of  nodes .  We  know  of  no  theorems  about  the 


existence  or  uniqueness  of  solutions  to  these  problems, 
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i  .  Closed  nonlinear  splines.  Now  consider  fitting  a  smooth 
closed  curve  through  a  set  of  prescribed  points.  We  will  express  this 
situation  by  utilizing  the  previous  development,  but  requiring  that  the 
points  Qq  and  Qn  be  coincident  at  an  n-th  prescribed  point,  and  that 
the  tangent  to  the  curve  be  continuously  turning  also  through  that  point. 
Thus,  for  some  integer  m  related  to  the  number  of  loops  in  the  curve. 


X0  =  Xn  '  y0  =  yn  '  0n  =  *0  +  ^  5 


(31) 


l  l 

nr  nr 

£  J  cos  6  ds  =  J  sin  0  ds  =  0  ; 


r=l 


r=l 


r-l 


’r-l 


)  (32) 


n  r 

J]  J  k  ds  =  2  m  tt 

r=1  l  . 
r-l 


The  deductions  from  (ig)  are  unchanged  from  those  described  heretofore, 

apart  from  the  contributions  at  s  =  l.  and  s  =  I  .To  obtain  a  local 

On 

minimum  of  the  integral  (14),  in  order  to  find  a  "smoothest"  closed  curve 
through  the  n  prescribed  points,  we  must  compare  curves  of  slightly 
different  total  arc  length,  and  this  can  be  achieved  by  selecting  the 
variations  and  to  be  unequal.  Since  the  tangent  to  the  curve 

prior  to  the  variation  is  continuously  turning,  and  that  after  the 
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variation  must  also  be,  the  variations  at  ln  and  t  must  satisfy 

On" 

6 e(tQ)  +  =  be(tn)  +  K(ln)btn  .  Co) 

It  is  not  correct  to  demand  that  be(l^)  =  °6(*n)  >  since  elements  of 
curved  arc  have  been  inserted  into  the  loop  in  superposing  the  variation. 
Since  o/Q  ,  Sf^  ,  56(*q)  and  6fl(in)  are  no  longer  independent, 

the  terms  arising  from  these  variations  in  (19)  must  be  combined  with 
(33)  to  deduce  the  natural  boundary  conditions  at  the  support  =  Qn  * 
At  the  boundaries  Iq  and  I n  ,  (19)  and  (22a)  give: 

2[/c(ln)5e(in)  -K(l0)be(t0)]  +  [K2(ln)  -2P(fn)]6ln 

-  [k2(Iq)  -  2P(I0)]6Iq  =  0  .  Ok) 

Eliminating  60(1^)  from  (33)  and  Ok)  gives: 

[k2(Iq)  +  2P(l0)]6f0  +  2l<(ln)  -  #c(f0)  ]&0(/n) 

+  [x2(ln)  -2 K(tQ)K(ln)  -PP(ln)]&ln  -  0  ,  (33) 

where  the  variations  bf^  ,  b6(i  )  and  of  can  now  be  considered 

0  *  '  n'  n 

arbitrary  and  independent.  Thus 

*(!0>  =  K(~lr)  Or>) 

and 

( X 2  +  2P)|  =  (K2  +  2P)  -  0  .  (37) 

0  n 

Thus  (29)  and  (30)  again  apply  throughout  the  spline.  Hence  the  natural 
boundary  conditions  for  (15)  yield  the  same  integration  constant  cr  =  0 
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in  (12)  for  the  closed  spline  as  for  the  free-ended  open  one.  However, 
for  the  closed  spline,  this  result  does  not  follow  from  the  least -constraint 
discussion  of  conservative  systems.  In  fact,  either  adding  or  removing 
an  element  of  arc  from  the  optimum  configuration  increases  the  strain 
energy  at  equilibrium  and  hence  exhibits  this  property  associated  with 
imposing  additional  constraint. 

7  .  Comments  and  examples.  When  curve  fitting  with  smooth  curves 
is  investigated,  the  variational  principle  (15),  utilizing  natural 
boundary  conditions,  calls  for  continuity  of  K  and  P  across  supports, 
as  well  as  the  continuity  of  6  prescribed  in  the  formulation  of  the 
problem.  Geometrical  discussions  of  the  problem  commonly  take  into 
consideration  only  continuity  of  6  and  K  ,  but  this  seeming  omission 
of  P  is  in  fact  automatically  taken  care  of  by  the  differential 
equation  (30),  satisfied  by  the  spline  in  each  span  between  supports, 
since  (30)  and  (29)  are  synonymous. 

The  variational  principle  (2)  will  yield  (23)  and  (10),  and  hence 
the  differential  equations  (12)  or  (11),  for  types  of  support  other  than 
the  least  constraining  one  treated  in  Section  4  above.  These  include, 
for  example,  pin  supports  which  prevent  sliding,  built-in  supports  which 
prevent  both  displacement  and  rotation,  and  a  fixed-angle  freely  displacing 
constraint.  In  general,  with  such  supports,  the  constants  c^  in  (23)  will 
not  be  zero,  and  will  change  from  span  to  span  along  the  spline,  so  that 
the  differential  equations  (IP)  or  (11)  govern  the  deflection  of  the 
spline  spans,  and  not  the  special  case  ( ;'0)  .  These  comprise  the  more 
general  elastica  curves  discussed  in  |Y|,  t’or  which  applied  forces  are 


i 


not  all  acting  in  the  direction  of  the  normal  to  the  spline  at  the  point 
of  application,  or  for  which,  in  the  closed  spline  case,  the  spline  does 
not  have  the  optimum  length  corresponding  to  (37)-  Note  that  in  the 
case  of  a  pin  support  must  be  zero,  and  when  rotation  is  prevented 

5fi(l  )  =  0  ,  and  it  is  such  conditions  which  modify  the  treatment  of  the 
previous  section. 

The  limiting  case  of  linear  splines  corresponds  to  beam  theory  when 
the  deflections  y  from  the  unstrained  spline,  considered  to  lie  along 
the  x-axis,  are  such  that  |dy/dx|  «  1  .  To  sufficient  accuracy,  x  can 

replace  arc  length  s  and  the  support  forces  can  be  considered  to  act  in 

the  y  direction,  and  then  the  longitudinal  force  P  is  zero  throughout. 
From  (10)  the  differential  equation  for  the  spline  then  takes  the  form 

d2x 

H  =  0  ,  (38) 

dx 

with  the  linear  approximation 

*  -  - 1  •  (59) 

dx 


This  immediately  leads  to  piecewise  cubic  polynomials  for  y  as  a  function 
of  x  .  The  variational  principle  for  linear  splines  is  that  they  minimize  (l) . 

Schweikert  [  11]  has  treated  linear  splines  under  tension,  in  which 
end  supports  supply  a  positive  longitudinal  force  P  ,  which  is  constant 
throughout  the  beam,  for  freely  sliding  constraints.  By  linearization 
of  (10)  it  follows  that 
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so  the  solution  between  successive  supports  takes  the  form 


y  =  c^  +  c2x  +  c^  cosh(ax)  +  c^  sinh(ax)  , 

where  a  =  V?  .  One  reason  for  introducing  tension  is  to  remove  extraneous 

points  of  inflection  of  the  interpolating  spline  curve.  The  variational 

principle  for  linear  splines  under  a  given  tension  P  is  that  they 

minimise  the  tota1  energy  of  the  system,  which  leads  to  minimizing 
x 

n  0  0 

J  [f"(x)  +  Pf*(x)  ]  dx 

*0 

among  all  functions  f  that  satisfy  the  constraints  and  have  continuous 
second  derivatives.  One  could  also  study  nonlinear  open  splines  under 
tension. 

The  theory  presented  heretofore  leads  to  seme  interesting  character¬ 
istics  for  particular  situations.  For  example,  both  for  the  open  spline 
with  minimum  constraints  depicted  in  Figure  2,  and  the  closed  spline  of 
optimum  length,  (29)  requires  that  P  be  zero  or  negative,  and  zero  only 
where  the  spline  arc  is  straight.  Thus,  whatever  the  geometry  of  the 
curve  being  fitted,  tensile  resultant  longitudinal  forces  will  never  occur 
(unless  they  are  imposed  at  the  ends) . 

Consider  now  fitting  a  closed  spline  through  the  vertices  of  an 
equilateral  triangle.  If  the  spline  is  bent  into  a  circle,  we  see 
from  (J)  that  M  is  constant,  whence  from  (5)  S  =  0  ,  and  from  (6) 

P  i  O  .  Hence  (29)  is  violated  by  a  circle.  To  satisfy  (29)  some 
additional  arc  length  must  be  added  to  produce  a  compressive  force  P  . 

The  '‘optimum"  spline  will  take  the  form  illustrated  in  Figure  4.  A  quali¬ 
tative  understanding  of  this  deduction  can  be  achieved  by  noting  that 
increasing  the  arc  length  for  a  given  angle  of  bend  tends  to  reduce  the 
contribution  to  the  integral  (l4),  just  as  adding  large  loops  to  a  spline 
configuration  permits  the  integral  (14)  to  be  reduced  towards  zero,  as 
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2 

mentioned  in  [3].  With  radius  R  ,  the  l/R  of  the  integrand  dominates 
the  2  rr  R  of  the  total  arc  length,  for  increasing  R  •  However,  for  a 
fixed  arc  length  and  total  angle  of  bend  (J  k  ds)  ,  the  contribution  to 
(lU)  is  a  minimum  when  K  is  constant.  Increasing  the  arc  length  of 
the  spline  in  Figure  4  from  the  circle  configuration  causes  a  variation  in 
curvature  which  tends  to  increase  the  integral,  offsetting  the  reduction 
associated  with  increase  in  arc  length.  The  latter  dominates  initially, 
to  yield  an  optimum  fit  illustrated  in  Figure  4. 

This  example  permits  an  assessment  of  the  interpolation  strategy 
expressed  in  (2),  since  one  might  regard  the  circumscribing  circle  as 
providing  a  more  natural  fit  through  the  vertices.  The  advantage  of 
increasing  the  arc  length  in  reducing  the  integral  (l4)  is  the  feature 
which  leads  away  from  the  constant -curvature  circle.  Inhibition  of  such 
a  tendency  can  be  achieved  by  imposing  a  penalty  on  increase  in  arc  length, 

for  example,  by  replacing  (2)  by 

l 

n  o 

5  J  (x  +  k)ds  =  0  .  (4o) 


Equation  (37)*  and  hence  (29)*  must  then  be  replaced  by 


} 


(41) 


so  that  for  this  simple  case,  choosing 


k  = 


'0 


* 


(42) 


where  is  the  curvature  of  the  circumscribing  circle,  yields  that 

circle  as  the  optimum  fit  according  to  (4o) .  Whether  such  an  approach 
could  be  generalized  is  an  open  question. 
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If  for  a  closed  spline  loop  passing  through  prescribed  points,  the 
arc  length  is  slightly  shorter  or  longer  than  the  optimum  length  given 


by  (29),  the  integral  (14)  will  be  larger  than  for  the  optimum  case. 

For  each  of  these  problems,  with  fixed  arc  length,  (2)  is  satisfied  by 
the  curve  form  assumed  by  the  spline.  An  illustrative  example  is  given 
in  Figure  5-  For  the  shorter  spline  loop 

\  +  p  >  0  (43) 


and  for  the  longer  one 

<  o  • 


(44) 


These  conditions  will  change  the  constant  in  the  governing 

differential  equation  (12),  which  will  apply  throughout  the  spline  with 
constant  c  ^  if  the  supports  are  freely  sliding  and  rotating. 

This  paper  has  treated  the  global  problem  of  spline  geometry.  The 
computation  of  spline  functions  to  approximate  the  spline  configurations 
considered  here  has  not  been  discussed  in  this  paper,  and  constitutes  a 
challenging  problem  in  numerical  analysis.  For  the  open  spline,  the 
curvature  at  the  first  support  is  zero,  so  that  only  the  angle  need  be 
determined  if  an  initial-value  approach  (the  so-called  "shooting  method") 
is  used  for  integration  of  the  spline  differential  equation  problem.  In 
the  general  closed  spline  case,  both  angle  and  curvature  at  a  support 
must  i  e  selected  for  an  initial-value  approach,  thus  posing  a  more 
cumbersome  problem.  For  the  problem  of  the  equilateral  triangle,  symmetry 
can  be  used  to  reduce  the  complexity  of  the  general  case.  However,  the 
work  of  Woodford  [12]  makes  it  seem  unlikely  that  shooting  is  a  good  way 
to  compute  splines. 
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Figure  2.  Spline  passing  through  supports. 


Figure  Forces  and  moments  on  a  spline  arc. 
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Figure  4.  Spline  fitted  through  the  vertices  of  an  equilateral  triangle. 
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